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Abstract 

We report the observation of Raman solitons on numerical simulations 
of transient stimulated Raman scattering (TSRS) with small group veloc- 
ity dispersion. The theory proceeds with the inverse scattering transform 
(1ST) for initial-boundary value problems and it is shown that the explicit 
theoretical solution obtained by 1ST for a semi-infinite medium fits strik- 
ingly well the numerical solution for a finite medium. We understand this 
from the rapid decrease of the medium dynamical variable (the potential 
of the scattering theory). The spectral transform reflection coefficient can 
be computed directly from the values of the input and output fields and 
this allows to see the generation of the Raman solitons from the numer- 
ical solution. We confirm the presence of these nonlinear modes in the 
medium dynamical variable by the use of a discrete spectral analysis. 



1 Introduction 

Stimulated Raman scattering (SRS) is a 3-wave interaction process with ex- 
tremely wide application in physics, especially in nonlinear optics This 
essentially nonlinear phenomenon couples two electromagnetic waves (pump and 
Stokes waves) to a two-level medium and is described by a simple system of par- 
tial differential equations j|. 

Such a system applies to different physical situations, depending on a phc- 
nomenological damping factor chosen to match the observed Raman linewidth. 
For instance, the steady state regime occurs when one neglects dynamical effects 
on the medium. Then the system becomes explicitly solvable in terms of the in- 
tensities and the result fits well the situation of strong damping and long pulses, 
such as in fiber guides When instead dynamical effects and damping are 
of same order, very interesting phase effects have been discovered in These 
have been interpreted as the manifestation of solitons which, in the spectral 
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transform scheme, would be related to discrete eigenvalues. Instead they are 
related to the continuous spectrum (they are not solitons) and named Raman 
spikes ||. Therefore the question of the observation of the Raman soliton is 
open since the original work || in 1975. 

When the damping term is much smaller than the dynamical response of the 
medium, and can be neglected the regime is called hyper-transient and applies 
to short duration pulses. Here the system possesses a Lax pair [|, [?], ||] and 
can be treated by means of the inverse scattering transform (1ST) generalized 
to boundary- value problems ||, |l0|]. The observation of the Raman soliton is 
an interesting open problem, especially since the boundary-value problem for 
TSRS has been completely solved on the finite interval in [fill with the essential 



result that, as expected from the works |12|, 13, 14 , the field q(x,t) universally 
evolves toward the self-similar solution. 

The commonly used TSRS model is obtained as a 3- wave interaction process 
where the group velocities of the 3 fields involved are considered equal. As 



shown in 15 , this assumption is an asymptotic limit and consequently the 
group velocity dispersion (GVD) has to be taken into account (through the 
spectral extension of the input laser fields), leading to a modified SRS system 



(cq. (2.2) below). The main consequence of this fact is that Raman solitons 
(discrete eigenvalues in the nonlinear Fourier spectrum) are effectively generated 
by SRS in the medium |||. 

The aim of this paper is to show that Raman solitons are indeed created 
in a finite length medium and that their generation is accompanied with poles 
in the nonlinear Fourier spectrum (spectral transform). We follow the model 
derived and solved in |ll| and establish the following results. 

1 - Based on the assumption of a semi-infinite medium, 1ST gives an explicit 
solution^] which will be shown to be an excellent approximation of the finite 
medium case. We explain this from the fast decay of the potential q{x,t) for 
large x. Such is not the case in the zero-GVD case for which q(x,t) decreases 
as x -3 / 4 typical of the self-similar solution jl2[ [Tl| . 

2 - The reflection coefficient p of the spectral transform (or nonlinear Fourier 
spectrum) is expressed in terms of the input and output field envelopes, allowing 
us to check the appearance of Raman solitons (the real valued single poles of p) 
on the numerical solutions. This also provides a means to observe the generation 
of Raman solitons in experimental data. 

3 - A recursion formula for computing the spectral transform of a finite set of 
data has been recently proposed Jl6[ . We have implemented this discrete spectral 
transform and confirmed the creation of these Raman solitons in a simulation 
on a finite length. A remarkable feature here is that the numerical spectral 
transform can be quite easily implemented, its computation is faster, and it 
provides more information than a conventional Fourier transform. 



lr The solution is not explicit for the finite length case: both results of |15| and of give 
the answer by solving a system of Cauchy-Green integral equations. 
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After presenting the model in section 2, we report the numerical observation 
of the Raman solitons in section 3. Section 4 explains this by an analysis of the 
medium dynamical variable. 



2 The model 

We present here the extension of the SRS model of || to the case where the 
dispersion of the group velocity is not neglected. We consider the classical model 
as in j|, well adapted to molecular Raman scattering, and for which the medium 
is schematically represented by a collection X of harmonic oscillators coupled 
to the electric field E through the polarizability of the medium. 



2.1 Basic equations 

When the polarizability depends on the frequency of the applied field, the group 
velocity of the electromagnetic waves becomes frequency dependent. While a 
very small group velocity dispersion (GVD) has no consequence on the evolution 
of the amplitude, it has nontrivial effects on the phase dynamics. In this context, 
it has been proved in []l5f , using multiscale analysis, that for the electric field 

E(*,0 =^"-""'> j dk a{k, x ,t) ,-<- 

and medium dynamical variable q(x, t)e l ^ Kx ~ n< ' + cc. where ui — L02 = ^ and 
fci — hi — K, the resulting model of transient SRS can be written in the retarded 
time t = t' — x/v 

d x a = qbe 2lkx , d x b = -qae- 2lkx , 

d t q = -g J dk abe- 2lkx , (2.2) 

where the overbar stands everywhere for the complex conjugate. Note the con- 
servation (x-independence) of the flux density \a\ 2 + \b\ 2 . 

The above model differs from the usual SRS system, corresponding to the 
zero GVD case 

d x a = q b , d x b = -q a a , d t q = -gab , (2.3) 

by the presence of the integral over all possible realizations of the phase mis- 
match k. 
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2.2 Boundary value problem 



The question of interest is the time evolution of a couple (a, b) of pump and 
Stokes pulses of time duration T sent into a medium of length £. Hence the 
domain of integration is 

x€[0,£\, te [0,1*1. (2.4) 

and it is necessary to prescribe the values of the fields on the two boundaries 
t = and x = 0. The medium is initially it = 0) at rest (all molecules in the 
fundamental state), hence the initial datum for the field q is 

q(x,0)=0, (2.5) 

independently of the reference frame []. 

The input [x = 0) light pulses are arbitrary functions of time t with a given 
spectral distribution around k = 0, namely 

a(k, 0, t) = A(k, t) , b(k, 0, t) = B(k, t) . (2.6) 

We shall be working here with the representative example treated in Jl5| which 
consists in assuming a common spectral lineshape for A and B as a normalized 
Lorentzian lineshape for the intensities. More precisely we set 

\A(k,t)f = \A (t)f -tj^-j , \B(k,t)f = \B (t)\ 2 -Ti^ , (2-7) 

7T k z + K z IT k -\- K 

where the input pulse profiles A (t) and B (t) are arbitrary. 

It is useful also to understand the scale invariance of the SRS system, in- 
cluding its boundary values. By a simple change of variables it can be easily 
shown that the system (2.2) (characterized by the parameters g and £) together 
with the input boundary values ( |2.6|) (|2.7[ ) (characterized by the parameter k) 
bears the scale invariance 

I — * at , g — > g/a , k — » n/a . (2-8) 

This invariance has been checked on numerical simulations as one way to test 
the code accuracy, and it works so well that we have obtained indiscernible 
pictures. 



3 Laser fields 



We describe hereafter the 1ST solution of the boundary-value problem (2.6) 



for the system (2.2) using the results of |15[, and compare it with the direct 



numerical solution in the finite length case. 



2 Note that a prepared medium would correspond to a state q given at physical time zero, 
that is on the characteristic t = —x/v. 
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3.1 Sketch of the direct problem 



The direct problem consists in defining the spectral transform from the data 
of the potential q(x,t) and the boundary values A(k,t) and B(k,t) (as time t 
appears everywhere as an external parameter, we shall forget it here). This is 
done by defining the following Jost solutions (for k € R) 



<p(k, x ) = i- / du(o#*.o , 

Jo 

cb(k,x) = [* d£q(£)<p(k,0 e 2lk(x -V , (3.1) 
Jo 

which are both entire functions of k vanishing as k — > oo in the upper half-plane. 

s 

These two functions then allow to define the reflection coefficient p(k) and 
the transmission coefficient r(k) by taking the limit x — > oo of ip and <f>e 2lkx for 
fcgK, namely 

■i POO p oo 

- = 1 - / dt; q(£)<Kk,Q , - = - / d£ q(0nk,0e 2lki ■ (3.2) 
T Jo T Jo 

The coefficient 1 Jr is clearly an entire function of k and one can show |H| that 
the reflection coefficient p(k) is meromorphic in the upper half-plane with a 
finite number of single poles related to the solitonic part of the solution q. We 
have also the following unitarity relation for keK 

H 2 -l + M 2 ■ (3.3) 

It is easy to prove finally that the vectors (ip, —<fie~ 2tkx ) and ((j>e 2lkx , Tp) solve 
the same differential equation as the vector ( a, b ) in (|2.2|). Then by comparing 
their values in x = we readily obtain from (|2j 



a = Aip + B(j)e 2lkx , b = B<p - A<\>e- 2lkx . (3.4) 
3.2 Output pump pulse 



From (p.2|), (|3.3| ) and (|3.4| ), the output \a(k,£, t)| is explicitly given for I — > oo 



by the expression 



|a(fc,oo,i)| 2 = — L-jj |A-pS| 2 , (3.5) 

i + \pY 



where A and B are the input data defined in (2.6) 



3 The relation with notations of is: <fi(k) = ip^^k) = tp 2 2(k), <j>(k) = — = < / 3 i2(' c )- 
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The main result of |15| is that the function p(k, t) is obtained by solving the 
following Riccati time evolution 

Pt = -p 2 C+[m*} - 2p C+[cj>] - C+[m] , p(k, 0) = , (3.6) 

where the functions m(k,t) and 4>(k,t) are given from the input data by 

7 7T 7 7T 

m = - g AB*, ^=- g {\A\ 2 -\B\ 2 ) , (3.7) 

and where denote the following Cauchy integral 

1 f + °° dC 
C t\f\ = -\ , (u.^ KO- (3-8) 



Consequently, for given inputs A(k, t) and B(k, t), the expression (3.5) gives the 



explicit asymptotic output pump intensity from the solution of the evolution 



(3.6). It is worth remarking that, at any given time, p possesses possibly a 
finite number of simple poles whose time evolution is given by the nonlinearity 
of the Riccati equation. 

For practical purpose, we choose in ( |2.6| ) the Stokes wave input seed as a 
portion e~ 7 of the pump wave, namely 

B (t) = A (t) e-i . (3.9) 

In that case the evolution (|3.6|) can be explicitly solved |L5| 



p{k,t) = 7JT7-, r, 3.10 

cosh(d(fc, t) - 7) 



6(k,t) = ^-, T(t) = -g{l + e- 2 r ) f d T \A (r)\ 2 . (3.11) 

K -f- IK 4 ,/q 



3.3 Numerical solution of finite length TSRS 

Our purpose is to understand how the above IST-solution can be used to model 



the solution of the SRS system (2.2) on a finite length. This can be done first 



in a qualitative way by writing the solution of (2.2) with boundary values ( |2.6| 
as the equivalent integral form 



a(k,x,t) \_f A(k,t) \ f x ( q(t,t)b(k,U)e 2zk Z 1 ( ~, 9 , 
b(k,x,t) J { B(k,t) J + J -q(t,t)a(k,£,t)e-™i I • 



The output a(k,£,t) will not differ much from its asymptotic value a(k, 00, t) 
if q(x,t) is sufficiently small for x > £. Consequently the behavior of q(x, t) at 



large x, is essential for the adequation of the formula (3.5) to real situations 



We will see in the next section that it is also crucial for the applicability of 
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the spectral method. We now proceed to show that the numerical solution of 



the TSRS system (2.2) on a finite interval is in excellent agreement with the 



theoretical expression for the output (3.5) 



We discrctizcd ( |2.2[ ) in both x and t using an order 2 Runge-Kutta method 
and advance via the following algorithm: 

1 - given a(k, x, t) and b{k 1 x, t) for all real k, compute q{x 1 1) by integrating the 
time-evolution of q at x for the initial datum q{x, t = 0) = 0, where the integral 
is calculated using the trapezoidal rule, 

2 - advance to a(k,x + dx,t) and b(k,x + dx,t), for all k, by integrating the 
differential equation for a and b, 

3 - go to step 1 with x — x + dx. 

The scheme is started at step 1 for x = 0. 

The quality of the computation is monitored by evaluating the relative error 
in the total flux \a(k,x,t)\ 2 + \b(k,x,t)\ 2 which is conserved by (2.2). In all the 
runs that are presented it remained smaller than 10 -5 . 

We chose as parameters 

^ = 80, .g = 0.5, 7 = 5, k = 0.2, (3.13) 
and the input pump pulse envelope is the Gaussian 



Ao(t) = cxp 



i-50 x L ' 



30 



(3.14) 



For this choice we found that dt < T/1000 and dx < re/2 gave stable results. 
Another point is that because of the Lorentzian line width we had to take a k 
interval of width w 80k in order to ensure that the integrals were normalized. 
To describe the strong oscillations present for the zero GVD system ( |2.3| ) we 
had to chose dt < T/2000 and dx < £/ 10000. 

A typical run with number of grid points in x, t and k (n x = n t = 1000; = 
500) takes about 2 hours CPU monoprocessor on a RS10000. We used the 
parallelism of the problem, i.e. the fact that the marching in x (resp. t) can be 
made in parallel for the loops in t and k (resp. x) and implemented the code on 
a Silicon Graphics SGI 10000 using the OpenMP software. This enabled a gain 
of a factor 8 or 10 in computing time depending on the number of processors 
used. 

Figure |l| show the pump intensity input |a(fc,0,i)| 2 and output \a(k,£, t)\ 2 
computed from the 1ST expression ( |3.5[ ) together with the numerical solution of 
the system ( |2.2[ ) for a length £ = 80 and four different values of the parameter 
k, in excellent agreement (a tiny discrepancy is only seen in the last picture). 

It is worth remarking that one cannot pursue a given numerical experiment 
to a longer length arbitrarily. We have stopped for instance the above calculation 
at length £ — 80 where the potential < 10~ 5 . Continuing the run to 

longer lengths would yield an increase of the potential which would then start to 
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oscillate. This is probably due to the system itself which amplifies the numerical 
errors in a drastic way at long lengths (note that from the scale invariance ( |2.8D 
longer length means larger Raman amplification). 



3.4 Raman solitons 



The function p of ( |3.1C ) has an essential singularity in k = —in and a set of 



single poles k n evolving in time, given by 

knit) = -IK + ^ , (3.15) 

for n G N. At time zero no pole is present and, as t evolves, poles move 
upward from — in and eventually reach the real axis at the times t n defined by 
Im(/c„(i„)) = 0, i.e. by the implicit expression 

r(*n) = - [ 7 2 + [(n+i)^] 2 ] . (3.16) 
7 2 

Note that t n = t_ n _i, hence the poles cross the real axis by pair. The corre- 
sponding positions on the real axis are then given by 

K 1 

k n (tn) = —k- n -l(tn) = Cn , Cn = ~(n + -) n ■ ( 3 - 17 ) 

7 2 

We plot in Figure the imaginary part of the poles k n (t) as a function of time 



for the parameter values (3.13) and n = 0, 1, 2 and 3. We observe that the first 
three poles cross the real axis at the positions Co = 6.310" 2 , £i = 0.19, C2 = 0.31 
and times to = 39.2, t% — 46.3 and t 2 = 59.6. The pole ks starts to evolve with t 
but cannot cross the real axis because of the finite duration of the pump pulse. 

As soon as these poles move to the upper half-plane, they generate a soliton 
component in the "potential" q(x, t). We will prove in the next section that this 
potential (the medium dynamical variable) is a continuous function of t when a 
pole crosses the real axis. 



3.5 The spectral transform from the output laser pulses 



The expressions (|3.4|) can be inverted to get the following expressions of the 
Jost solutions in terms of the (physical) fields a{k, x, t) and b{k, x, t) and of the 
input pulses A(k,t) and B(k,t) 



aA + bB 



\Af 



B 



0e 



2 //,•./• 



aB - bA 



\Af 



B 



Then the definitions (3.2) lead to the following formula 

bA - aB 



p{k,t) 



aA + bB 



(3.18) 



(3.19) 
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which gives the spectral transform p(k, t) in terms of the output pump and 
Stokes fields. This function for k real must become singular at the two points 
±£ n each time a soliton k n is created and it is used now to prove the generation 
of Raman solitons in numerical experiments. 



For the input pulses given in (2.7) and the particular choice ( |3.9| ), we readily 
obtain from the above 



ae 



(3.20) 



We use this expression to estimate pi from the numerical solution (a, b) at 
x — £ < +00 and compare it to the p obtained from the 1ST (3.1C) for the 
parameters (3.13). Notice that \p(— k)\ = \p(k)\ so that we will only present 
positive values of k. 

Figure || shows the function \p\ of (3.1C) in full line and the numerical result 
obtained from (3.2C) in x = I the in dashed line, for the 3 values k = Co, £i> C2 
and 0$. Both expressions are very close and as expected p(( n ) is singular for 
t = t„ for n < 2, while it is regular for n = 3. In the first picture of figure 
||, the dashed line actually represents the value at k = instead of k = £q. 
This is the only noticable discrepancy between analytic asymptotic formula and 
numerical/experimental expression, resulting from the finiteness in x of the data 
q(x) (the wave length A = 27r/Co ~ 104 > I). 



4 Medium 

In the previous section we have seen that the finite length solution is in good 
agreement with the asymptotic behavior ( |3.5| ) given by the 1ST on the semi- 
infinite line. We now proceed to justify this fact by analyzing the behavior of 
the medium dynamical variable q(x,t). 



4.1 Sketch of the inverse problem 

1ST furnishes the solution q(x,t), called medium dynamical variable, by the 
expression [[15) 

q(x,t) = 2iiP w (x,t) (4.1) 

where V^ 1 "* is the coefficient of 1/fc in the Laurent expansion of the function 
ip(k,X,t) solution of the following Cauchy-Green coupled system 

{k) = 1 + J_ f + ° dx e 2iX X _ y pjknMkn) e 2f knX 

2i7r J_ 00 A — k — iO k — k n 

2«7i" J -00 X — k + iu k — k n 
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Here p n {t) are the N residues of p(k,t) at the poles k n {t) that are in the upper 
half fc-plane (if any). Note that the solution cp corresponds to the solution of 
([□]) and that if;(k) =~<t>(k). 

Apparently the generation of a soliton is followed by the adjunction of a term 
in the equation for tp, and consequently also in the expression of q, leading to 
a possible discontinuity for q. We will now proceed to show that this is not the 
case. This situation is particular to the spectral problem on the semi-line where 
the continuous spectrum (value of p(k) on the real axis) is not separable from 
the discrete spectrum (poles of p in the upper half complex plane) . This is due 
to the motion of the poles of p that can cross the real axis. On the contrary, 
in the full line case, solitons (discrete spectrum) can exist without radiation 
(continuous spectrum). 



4.2 Continuity 

From now on when the spectral parameter k is generically complex we shall be 
using boldface letter k. From what precedes, at t — to a pole crosses the real 
axis at fc = fco from the lower half plane. We consider here the limits t — > to ± 
and we write for k S M 

Pdk)=p(k,t ±0) = ] ^- (4.3) 

where i?o(k) is analytic in a neighborhood of k . Note that 

p + (Jt)=p-(Ji) + 2mp 6{k-k ) (4.4) 
where po = Ro(ko) is the residue of p at the pole k — fco. 



Assuming for simplicity no other pole, we have from (4.2) at t = to — 

m 1,40-0) = -?- r dX MAMA Mo -0) £ _ 2lXx 
2«7T j-oo A — k + iv 

which is well defined at k = fco since the distribution (A— ko+i0)~ 2 is meaningful. 



On the contrary at t = to+0, expression (4.5) diverges due to the distribution 
(A — fco + iO) _1 (A — fco — i0) _1 . Since now there is a pole in the upper half complex 
plane, we should take the limit Im(k) — > in the complete expression in (4.2). 
i.e. 



1 / >+0 ° 
t/j(k,x,t o + 0) =— / dX 
2«tt J_ OQ 



p + (X)<p(X,x,t o +0) 



A-k 



e 



+ /w(fcp,:Mo + 0) c 2 %kax ( 4 
k - fc 

for Im(k) < 0. We obtain 

^(Mo-0) =^(Mo + 0). (4-7) 
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and consequently from (4.1) 

g(*o-o) = g(t + o). 



(4.8) 



4.3 Asymptotic behavior 

We consider here for simplicity the case when solitons are not yet present. Then 



from (4.1) and (4.2), q can be reconstructed in terms of the reflection coefficient 



p and the Jost solution ip by means of 

l(x) = / dk p(k)<p(k,x)e- 2tkx . (4.9) 

This formula yields directly an information on the behavior of q at large x. 

In fact, if p(k)(p(k, x), . . . , {p(k)(p(k, x)} are continuous and tend to 

for k — > oo, and if {p(k)(p(k, x)} £ L(WL), we obtain by repeated integration 
by parts 

q(x) = (J^) J dk d^{p(kMk,x)}e-^. (4.10) 

Hence finally 

x — ► oo x n q(x) -» . (4.11) 

As mentioned in the introduction, such an asymptotic behavior is not found 
for the potential qo(x,t) which would be obtained from TSRS with zero GVD 
(p75|). Indeed, in that case the medium initially at rest evolves universally to- 
wards the self-similar solution |ll[] which behaves as a; -3 / 4 as found in fll2|] . 
More precisely we have 

Q~s(t)C 3/i [acos(he /2 +l3) + 0(C 1/2 )] , £ = s(t)x, (4.12) 

where s(t) is a function of t determined by the boundary conditions and where 
a, and h are arbitrary constants. 

Figure U presents the decay of |g| and |go| respectively in dashed line and full 
l ine a s a function of x for t = T/2 = 50 in a linear-log plot and the parameters 
( |3.13 ). The exponential decay of \q\ can be seen for all values of t and follows 



approximately e~ KX as expected from the ana lysis. This observation justifies 
the fact that the finite length evolution of ( |2.2| ) is very close to the asymptotic 
expression (^) given by 1ST. 

On the contrary the slow decay of |<7o| as a; -3 / 4 makes it impossible for the 
integrals of (|]l]) to exist so that the scattering theory must be revisited in this 
case. 
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4.4 Numerical spectral transform and Raman solitons 



By numerical integration of the SRS system (2.2) with a spatial grid of dimension 
h, we get, at each value of time, a set of L discrete data q(n,t) (with x = nh 
and I — Lh), that we now analyze by means of the spectral transform. This is 
easily done, using the results of |l6| which, adapted to our notations, read 

p(k,t) = Ro(C,t) , (4.13) 

where Rq is obtained from the following inverse recursion 

r , , R m - hq{m-l,t) C m ~ l 

R L = C L h q (L,t), ^-^ ^V^m-MK-^ ' (414) 

and where the parameter £ is related to k by 

C = e 2lkh . (4.15) 

It is then a simple task to use this recursion relation to compute the numer- 
ical spectral transform i?o(CiO an d compare it to the asymptotic theoretical 
expression ( 3.10 ) of p(k,t). In Figure [| we plot \p(k,t)\ where one can clearly 



see the three singularities (Cn, t n ) for n — 0, 1, 2 corresponding to the three poles 
of p. 

The non linear spectral transform is easier to compute, faster and yields 
more information than the standard Fourier transform which in particular does 
not present any singularity (as we checked numerically). 



5 Conclusion 

This work demonstrates that the 1ST solution of transient SRS on the semi- 
infinite line furnishes a very accurate model for the solution on a finite domain. 
The accuracy stands not only for the output laser pulses intensities but also 
for their phases as shown in a spectacular way by the generation of the Raman 
solitons. It is shown that the question of the experimental observation of Ra- 
man solitons is solved by a convenient combination of the output (and input) 
laser profiles. We also performed a numerical nonlinear spectral analysis which 
resulted to be not only rich of information but also quite easy to implement and 
faster than the usual FFT procedure. 

Finally we mention that it is difficult to compare these results with those 
obtained in the zero GVD case ( |2.3| ). The problem is that in the limit k — > 0, 
the spectral transform gets an essential singularity in k = —iO which sends poles 
to the upper half-plane. To discuss the number and time location of such poles 
would require not only a careful study of the singular limit k — > 0, but also to 
reformulate 1ST on the half-line as we have seen that the potential q(x) does 
not decrease fast enough as x — > oo. 
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Figure 1: Plot, as functions of time, of the output pump intensity \a(x = 
£,k,t)\ 2 from the theoretical expression (3.5) (full line) and the numerical so- 
lution of (2.2) (dashed line) for four different values of k. The large gaussian 
curves (dashed line) represent the input pump intensities \a(x = 0,fc,i)| 2 . 
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Figure 3: Time evolution of the reflection coefficient |/c(Cn)| for n = 0, 1, 2 and 
3 obtained from the inverse scattering theory fl3.1Q ) (full line) and from the 
numerical solution of (|2.2|) using formula (3.20) (dashed line). 
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Figure 4: Evolution of \q\(x, t — 50) and l^oK^, t — 50) in linear-log coordinates. 
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Figure 5: Evolution of |p|(fe,t) obtained from the discrete spectral transform 
(pi. 
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